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We report the derivation and implementation of orbital optimization algorithms for the active space decom¬ 
position (ASD) model, which are extensions of complete active space self-consistent field (CASSCF) and its 
occupation-restricted variants in the conventional multiconfiguration electronic-structure theory. Orbital rota¬ 
tions between active subspaces are included in the optimization, which allows us to unambiguously partition the 
active space into subspaces, enabling application of ASD to electron and exciton dynamics in covalently linked 
chromophores. One- and two-particle reduced density matrices, which are required for evaluation of orbital 
gradient and approximate Hessian elements, are computed from the intermediate tensors in the ASD energy 
evaluation. Numerical results on 4-(2-naphthylmethyl)-benzaldehyde and [3e]cyclophane and model Hamil¬ 
tonian analyses of triplet energy transfer processes in the Closs systems are presented. Furthermore model 
Hamiltonians for hole and electron transfer processes in a/iri-[2.2](l,4)pentacenophane are studied using an 
occupation-restricted variant. 


I. INTRODUCTION 

Electron and exciton dynamics between chromophores are 
ubiquitously found—from photosynthetic complexes in bio¬ 
logical organism^to organic semiconductors for solar energy 
conversion. 2 They are often modeled by quantum dynamics 
methods based on the quantum master equation, 3 which re¬ 
quires so-called model Hamiltonians in the diabatic represen¬ 
tation as an input. The model Hamiltonians are the reduced 
dimensional representation of the electronic degree of free¬ 
dom in such systems, consisting of the energies of the dia¬ 
batic states and the interaction strengths between these states. 
Although accurate model Hamiltonians are of essential im¬ 
portance in achieving predictive simulations of electron and 
exciton dynamics, there have been only a few studies to de¬ 
velop accurate electronic structure methods to compute them 
from first principles.® 1 

We have recently introduced the active space decompo¬ 
sition (ASD) method to provide accurate model Hamiltoni¬ 
ans. The ASD model takes advantage of molecular geome¬ 
tries to compress active-space wave functions of molecular 
dimers, 12 !® in which wave functions are parameterized as?- 2 

r?> = X (i) 

u 

where A and B label monomers and I and J label monomer 
states. The monomer wave functions (|<I) 4 ) and O?)) are 
determined by diagonalizing the respective monomer active 
space Hamiltonians (i.e., // 4 |<F 4 ) = Ef |0 4 )) in an orthogonal 
active subspace (hereafter referred to as an ASD subspace). 
The expansion in Eq. 0 is exact when each I and J en¬ 
tirely span the corresponding monomer space, and converges 
rapidly with respect to the number of monomer states included 
in the summation. 12 An extension to more than two active 
subspaces has also been reported by the authors. 15 Further¬ 
more, the dimer basis states have well-defined charge, spin, 
and spin-projection quantum numbers on each monomer, al¬ 
lowing us to extract model Hamiltonians for electron and ex¬ 
citon dynamics through diagonalization of diabatic subblocks 
of a dimer Hamiltonian matrix. 13 This methodology has been 


applied to singlet exciton fission dynamics in tetracene and 
pentacene crystals. 14 The resulting model Hamiltonians have 
been used to benchmark approximate methods in more recent 

studies.® 

However, our ASD method was not applicable to covalently 
linked chromophores in previous works. This was mainly due 
to the difficulty in finding appropriate ASD subspaces. In the 
original formulation of ASD, one must define localized ASD 
subspaces from localized Hartree-Fock orbitals.® Although 
well defined (and even automated) for molecular dimers and 
aggregates, this procedure gives rise to an ambiguity in the 
definition of the ASD subspaces in covalently linked chro¬ 
mophores owing to orbital mixing. This problem has ham¬ 
pered applications of ASD to these systems. 

In this study, we remove this ambiguity in the subspace 
preparation by deriving and implementing orbital optimiza¬ 
tion algorithms for the ASD model. The orbital optimiza¬ 
tion guarantees that (when converged) identical results are ob¬ 
tained regardless of the choice of initial active orbitals. As 
shown in the following, the orbital optimization procedure 
naturally leads to localized ASD subspaces without deteri¬ 
orating the diabatic structure of the dimer basis function in 
the original ASD method. The preservation of locality of the 
ASD subspaces during orbital optimization is analogous to the 
fact that, for a fixed number of renormalized states, ab initio 
density matrix renormalization group algorithms using local¬ 
ized orbitals provide the lowest energies. 17 A quasi-second- 
order algorithm is employed,® which is akin to conventional 
complete active space self-consistent field (CASSCF) and its 
occupation restricted variants (RASSCF). 19 In the following 
our orbital-optimized models are referred to as ASD-CASSCF 
and ASD-RASSCF depending on the underlying monomer 
active-space wave functions. 

H. THEORY 

A. Orbital optimization model 

In ASD the energy of a dimer state is a function of molec¬ 
ular orbital (MO) coefficients (C), configuration interaction 



2 


(Cl) coefficients within each monomer (c^ with D labeling 
Slater determinants for X — A, B), and ASD coefficients [Uu 
in Eq. (QJ], i.e.. 


£asd = E(C,4,4„U u ). (2) 

We first form the monomer basis by solving configuration in¬ 
teraction within each ASD subspace 

H a |®?> = Ef |<D?>, (3) 

H b \® B j ) = E B \® B ), (4) 

where |®^> = Ed c o' | D) and so on. Using the monomer basis 
states, the dimer Hamiltonian matrix elements are computed 
as 


= (-1)* (5) 

(n 

which can be evaluated without constructing the product ba¬ 
sis functions explicitly. Here the indices £ and ;/ are opera¬ 
tors acting on monomer A and B, respectively, whose product 
corresponds to either one- or two-electron operator of the re¬ 
arranged Hamiltonian; h^ n is the molecular orbital integral; 
and, (—1)^ is a phase factor due to the rearrangement. The 
monomer intermediates T AJ! and F BJJ are defined by 

T? r = ( 6 ) 

r BJJ ' = (® B \E n |®*>. (7) 

Diagonalization of the dimer Hamiltonian matrix gives ASD 
energies and wave functions. The details of the algorithm are 
given elsewhere. 12 

We introduce a stationarity condition with respect to vari¬ 
ations of MO coefficients to define orbital-optimized ASD 
models: 


<5 £asd 

dK rs 


( 8 ) 


where r and s label any molecular orbitals, and k is the anti¬ 
symmetric matrix that parameterizes the MO coefficients. 


C = Cinit exp(tf). 


(9) 


The non-redundant elements of the rotation matrix k are de¬ 
picted in Fig. □ The two-step optimization procedure (i.e., 
alternating solutions of ASD and orbital updates until conver¬ 
gence) gives a unique energy that is close to the variationally 
optimal one. This quasi-variational behavior is explained as 
follows. The total derivatives of the ASD energy with respect 
to variations of MO coefficients are 


dE ASD _ SEasd v' 1 <3£asd dUu 

dK rs dK rs L-f dUiJ dKrs 


+ 2 


dEAsp dcp 1 
dc^ 1 dK rs 


z 


SEasd dc\ g/ 

dc B D f dKrs ’ 


( 10 ) 
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FIG. 1. Non-redundant elements of the rotation matrices k for 
CASSCF (left) and ASD-CASSCF (right). 


where respective normalization conditions are implicit. First, 
it follows from the ASD procedure that the stationary condi¬ 
tion with respect to Ujj, 


d^ASD 

dU u 


(11) 


is automatically satisfied. Second, the following equations are 
approximately satisfied, 


5£asd _ „ 5£asd 

3C 4 ~ dc B 

ul d uc d , 


( 12 ) 


The residuals of these derivatives are related to the differ¬ 
ences between the ASD energies computed by lowest-energy 
monomer states and variationally optimal ones (with the same 
number of states in the summation), which are expected to be 
much smaller than the energies themselves in the cases where 
ASD is a good approximation. As a result, one obtains that the 
total derivative [Eq. (|T0|] is approximately zero, making the 
orbitals almost variationally optimal. In the cases where con¬ 
vergence is not achieved under this approximation, we fix cl 
and c b d , after several iterations (but not Uu)\ when c^ and cf y 
are fixed, the last two terms in Eq. ( |T()| that introduce numer¬ 
ical noise become exactly zero because the derivatives of c A ) 
and c B y are zero, and hence, the minimum in this constrained 
frame can be found by standard minimization algorithms us¬ 
ing the partial derivative of energy with respect to orbital ro¬ 
tations [i.e., Eq. <01- After the convergence ASD-CI is per¬ 
formed once again using the optimized orbitals. The solutions 
are obtained by iterating this procedure until consistency be¬ 
tween orbitals and c A D and c B D , is achieved. 

Lastly, we note that the orbital-optimized ASD energy be¬ 
comes equivalent to the conventional CASSCF and RASSCF 
energy when all the charge and spin sectors are included in 
the ASD expansion. In Sec. [HI] we show that the energy con¬ 
verges rapidly with respect to the number of states in the ASD 
expansion. 


B. Orbital gradients and approximate Hessians 

In addition to orbital gradients with respect to rotations be¬ 
tween closed, active, and virtual subspaces for which gradi¬ 
ent elements are knownj^we consider rotations among ASD 
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subspaces A and B. The energy gradients between orbitals in 
active space A (unbarred) and those in active space B (barred) 
are 


SEasd 2(F - F- ) 

dKij { jl ' jh 

(13) 

Fy = y t r j k Fi k + y r jkMl (ik\lm), 

(14) 

k klm 


F- u = h- u + y [2(7;|oo) - (io\oj)]. 

(15) 


O 


where i, j, k, l, and m label active orbitals (the summation 
indices k, l, and m run over orbitals on both A and B), and o 
labels closed orbitals. I and Vi L ki are the standard spin-free 
one- and two-particle reduced density matrices, respectively. 

In the quasi-second-order optimization algorithm 18 approx¬ 
imate Hessian elements are required as an initial guess. We 
use the following formula for approximate diagonal Hessian 
elements for the inter-subspace active-active rotations, 


d 2 £ASD 

dxijdKij 


~ 2[rs(^ + Ffj) + FjjiFL + Fp 
- Fv-Fjj-TTyF /.], 


(16) 


where F4 = Yjki^ki[(ij\kl) - \(jl\kj)]- This formula is derived 
using similar approximations to those in Ref. [18] in which the 
two-electron integrals are approximated by the diagonal ele¬ 
ment of Fock-like one-electron operators and active orbitals 
is approximated to be completely filled for some terms. This 
procedure eliminates the needs for the computation of two- 
electron integrals with more than one active indices. 

One- and two-particle reduced density matrices in the ASD 
model, Tij and T,-^/, are computed from 1)- intermediate ten¬ 
sors as using Eqs. 0- 0, and {7) as follows, 

(17) 

//' 

ff = J] UuUrrT*' 11 ', ( 18 ) 

jj' 

which after index reordering give F, ; and T,,-*.;. Because T k 
intermediates are already computed in the ASD energy evalu¬ 
ation, additional costs for computing density matrices are neg¬ 
ligible. Note that density matrices have permutation symme¬ 
try; for instance. 


t ijjd - r 'jij k - r kl jj - r lk ji . (19) 

The density matrix elements with all indices belonging to the 
same monomer are computed separately for computational ef¬ 
ficiency as 


m - X X {^kXj^j), 

J per 

(20) 

i]M= x yimyirim. 

(21) 


I per 


using rotated monomer states |€>^) = 2/ Ujj\<t>^) and |®f) = 
Yj.i Uij\A> B j). The use of rotated monomer states reduces the 
cost of calculating these elements from quadratic to linear 
scaling with respect to the number of monomer states. 


C. State averaging and model Hamiltonians 


As in traditional CASSCF, state averaging is used to op¬ 
timize orbitals for systems that involve multiple electronic 
states, in which £ a sd in Eq. ([8]) is related by an averaged en¬ 
ergies of states of interest (labeled by K), 

^ASI) = ~ X ^ASD- ^2) 

K= 1 

This means that all the density matrices r,- ; - and T,^/ that ap¬ 
pear in the orbital gradient and approximate Hessian elements 
have to be replaced by state-averaged counterparts. 



K= 1 


(23) 

(24) 


The rest of the algorithm remains identical. 

As discussed in the introduction, ASD-CASSCF is devel¬ 
oped to realize calculation of diabatic model Hamiltonians 
for electron and exciton dynamics. The procedure to extract 
model Hamiltonians from ASD dimer Hamiltonian matrices 
has been detailed in Ref. [TJj In essence, we form a diabatic 
model state by diagonalizing each diagonal-subblock of the 
dimer Hamiltonian corresponding to appropriate charge and 
spin quantum numbers. The diabatic wave functions are writ¬ 
ten as 


\K A n B )= X (25) 

IeQ. A Jen B 

in which and fig are collections of monomer states with 
given charge and spin quantum numbers. The diagonal el¬ 
ements are the eigenvalues in the subblock diagonalization. 
The off-diagonal Hamiltonian matrix element of the diabatic 
states are defined as 


Hkk' = ). (26) 

When computing model Hamiltonians with orbital optimiza¬ 
tion, it is convenient to use the averaged energy of the diabatic 
states in the model Hamiltonian. This procedure gives iden¬ 
tical energies to the state-averaged calculations with the adia¬ 
batic states obtained by diagonalizing the model Hamiltonian, 
since the trace of a Hamiltonian matrix is invariant under the 
unitary transformation. 


HI. RESULTS AND DISCUSSION 

We applied ASD-CASSCF to the benzene dimer in its sand¬ 
wich configuration with a separation of 4.0 A. The cc-pVDZ 
basis sePil was used, and singlet (S), triplet (T), and charge- 
transfer (CT) monomer states are included. The CAS (12,12) 
active space for a dimer consists of 12 n electrons distributed 
in the n orbitals, which is decomposed to two 6-orbital ASD 




ASD-CASCI A S 

ASD-CASSCF -©- S 


4 


S+T S + T + CT 

S +T -^S + T + CT 



FIG. 2. Errors in (a) the ground-state energy and (b) the energy 
splitting between the ground and first excited singlet and triplet states 
of a benzene dimer as functions of the number of states in each charge 
and spin sector in the ASD expansion ( N ) computed by ASD-CASCI 
and ASD-CASSCF. 


subspaces. First, the ground state energies and the first 
few excitation energies computed by ASD-CASCI and ASD- 
CASSCF as a function of the number of states in each charge 
and spin sector in the ASD expansion ( N ) are presented in 
Fig-0 The errors are calculated relative to the conventional 
CASCI (-461.547119 E h ) and CASSCF (-461.584702 E h ) 
energies. ASD-CASSCF and ASD-CASCI converge to the 
conventional CASSCF and CASCI with full CAS(12,12) at 
large N. The Hartree-Fock orbitals were used in the CASCI 
calculations. For the ground state, the convergence of ASD- 
CASSCF is slightly slower than that with ASD-CASCI, al¬ 
though they are both almost exponential. The energies con¬ 
verged to 0.1 mZin with as small as /V = 8. We also note 
that ASD-CASSCF with a more diffuse aug-cc-pVDZ basis 
set gave almost identical errors for given N. See the support¬ 
ing information for details? 1 The excitation energies show 
negligible variance with respect to the number of monomer 
states. In all cases, the inclusion of CT states is essential to 
obtain numerically near-exact energies. The T-T contributions 
were negligible for the ground state, since the configuration 
in which the two monomers have anti-parallel triplets lies too 
high in energy to interact with the neutral S-S configuration. 
Orbital optimization increased and decreased the CT contri¬ 
butions to the ground and excited states, respectively. 

The convergence of orbital optimization in ASD-CASSCF 
is compared to that in CASSCF and is shown in Fig. [3] Three 



FIG. 3. Convergence of orbital optimization in ASD-CASSCF and 
the standard CASSCF for a benzene dimer (see text for computa¬ 
tional details). A quasi-second-order algorithm was used. 




N 

FIG. 4. Errors in the ASD-CASSCF ground-state energies as func¬ 
tions of the number of states in each charge and spin sector in the 
ASD expansion ( N ) for (a) 4-(2-naphthylmethyl)-benzaldehyde, and 
(b) Gelcyclophane. 


states were state-averaged in the calculations. We found that 
the convergence behaviors of ASD-CASSCF and CASSCF 
were nearly identical, corroborating that active-active rota¬ 
tions do not deteriorate the numerical stability of orbital op¬ 
timization. The total ASD-CASSCF calculation (7 iterations) 
took about 75 and 210 seconds with N = 8 and 32, while 
the standard CASSCF took about 3300 seconds using 16 CPU 
cores of Xeon E5-2650 2.00 GHz. 

Next, the energies of 4-(2-naphthylmethyl)-benzaldehyde 
(M) and [3f,]cyclophane (CP) are presented in Fig.[4]as func¬ 
tions of the number of states in the ASD decomposition 
(see Refs. j22. [23] and [24 for experimental studies on these 
molecules). The geometries were optimized by density func¬ 
tional theory using the B3LYP functional and the 6-31G* 
basis set.® The optimized structures of these molecules are 





























5 



TABLE I. Diabatic coupling \Hjf\ (in meV) and triplet energy trans¬ 
fer rate k (in s -1 ) for the Closs systems. The factors in the rate ex¬ 
pression other than the diabatic coupling were taken from Ref. [6] 


Molecule 


ASD-CASSCF 

CIS-Boyr 


fcexpf| 

\H, f \ 

k 

\H if \ 

k 


1,3-C-ee 

3 

0.863 

7.0 x 10 9 

1.5 

2.1 x 10 10 

7.7 x 10 9 

1,4-C-ee 

4 

0.363 

1.6 x 10 9 

0.56 

3.9 x 10 9 

1.3 x 10 9 

2,7-D-ee 

5 

0.076 

6.9 x 10 7 

0.17 

3.5 x 10 s 

9.1 x 10 7 

2,6-D-ee 

6 

0.010 

1.1 x 10 6 

0.020 

5.0 x 10 6 

3.1 x 10 6 


a Values based on configuration interaction singles diabatized by Boys 
localization, taken from Ref.|6] 
b Number of cr-bonds between the donoi^and acceptor. 
c Experimental values taken from Ref. ns 


transfer is given b>0 


FIG. 5. Structures of 4-(2-naphthylmethyl)-benzaldehyde (M), 
[36]cyclophane (CP), and aut/-[2.2](l,4)pentacenophane (PP). 


k = -j~\H if \ 2 
n 


1 


4 tiJcbT A 


exp 


(/! +AG°) 2 \ 

4 Ak B T )’ 


(27) 


shown in Fig. [ 5 ] The def2-SVP basis set 27 was used in 
the ASD calculations. In this example, ASD-CASCI can¬ 
not be used for these molecules; therefore, only the results 
with ASD-CASSCF are shown. The active space for M con¬ 
tains all the ^--electrons from benzaldehydyl and naphthyl 
groups to form CAS(18,18), which is decomposed to a prod¬ 
uct space of (8 ± q, 8) and (10 + q, 10) where q (= 0,1, ■ • •) 
is the number of electrons transferred in the dimer basis. The 
construction of ASD subspaces for CP is analogous to that 
for a benzene dimer above. For CP, the CASSCF energy 
(-1156.239888 £),) is used as a reference, while the ASD- 
CASSCF energy with N = 128 (-764.327532 E h ) is used for 
M since the exact dimer calculations with CASH8,18) calcu¬ 
lation are not feasible. The large 7r-overlap of CP due to con¬ 
strained geometry yields the slower convergence with respect 
to the number of states and the pronounced polarization con¬ 
tributions to the energies in comparison to M; however, with 
N = 32 energies for both molecules are converged to less 
than 0.1 mEh- The errors presented in Fig. [4] are dominated 
by that from the ASD approximation to the active Cl, and not 
by variations of orbitals with respect to N. To verify this, we 
computed the ground-state energy using ASD-CASCI with or¬ 
bitals obtained by N — 8 ASD-CASSCF. The convergence 
was found to be almost identical to that of ASD-CASSCF, at¬ 
testing that the change of orbitals is a minor contribution to 
the ASD errors. 

We then turn to the triplet energy transfer processes of the 
donor-bridge-acceptor systems pioneered by Closs and co¬ 
workers, 2223 which have been extensively studied theoreti¬ 
cally by Subotnik and co-workersP® The donor and accep¬ 
tors are benzaldehyde and naphthalene, respectively, and the 
bridge is either cyclohexane or frans-decalin, each of which is 
rigid and saturated. In Marcus theory, the rate of triplet energy 


where AG° is the driving force, A is the reorganization energy, 
and Hip is the diabatic coupling between the two diabatic 
states, defined in Eq. ( |26[ ). The triplet energy transfer is medi¬ 
ated by charge-transfer states, whose contributions are effec¬ 
tively treated using the quasi-degenerate second-order pertur¬ 
bation theory in our model. The diabatic coupling now reads 13 


Hip 


- Hip + 


\ z 

Z£ model 


/ HizHzf 
\E,-Ez 


HizHzf \ 
E f - E z ) 


(28) 


Here, the second term refers to the perturbative correction to 
the direct coupling, and the sum runs over all diabatic dimer 
states Z that are not included in the model spaces. We use d ge- 
ometries optimized by B3LYP with the 6-31G* basis set. 25 ® 
Two triplet states were included in the state averaging, and the 
def2-SVP basis set 27 was used in the ASD calculations. The 
active spaces are the same as that for M [i.e., decomposed 
CAS(18,18)], and 64 states were included in each charge and 
spin sector of the ASD expansion. Table [I] compiles the cal¬ 
culated diabatic coupling and rates from the state-averaged 
ASD-CASSCF orbitals. The mediated coupling values are 
tabulated; corresponding direct coupling values are 0.074, 
0.040, 0.0078, and 0.00033 in meV for 1,3-C-ee, 1,4-C-ee, 
2,7-D-ee, and 2,6-D-ee, respectively. 

Combining the diabatic coupling elements with the prefac¬ 
tors reported by Subotnik and co-workersP the triplet energy 
transfer rates were calculated using Eq. Our coupling 

values are roughly half of those obtained by Subotnik and co¬ 
workers based on the diabatization of configuration interac¬ 
tion singles (CIS) adiabatic states via Boys localization. The 
calculated rates are in good agreement with the experiments. 
For more consistent comparison to experimental results, one 
needs to determine the prefactors using ASD-CASSCF. It is 
also noted that we obtained the diabatic coupling of 0.54 meV 
for M, one of the Closs systems, in good agreement with 0.56 
meV by Subotnik, although the Condon approximation is not 
valid for this molecule owing to the flexible methylene bridge. 
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FIG. 6. Initial (mesh) and converged (filled) active semi-canonical orbitals for 1,3-C-ee from state-averaged ASD-CASSCF. 


Figure [6] visualizes the initial and the optimized active or¬ 
bitals for 1,3-C-ee to show that the diabatic nature of dimer 
product functions is retained during the orbital optimization 
procedure including active-active rotations [see Eq. m- 
The shapes of the n* orbitals slightly changed during the opti¬ 
mization, while the 7r-orbitals remained almost identical. The 
locality of the optimized orbitals should be ascribed to the fact 
that the ASD expansion [Eq. (|TJ] is most accurate with local¬ 
ized ASD subspaces. The delocalization of the active orbitals 
to the bridge site is suppressed in these examples since the 
active space is constructed with only 7r-orbitals and the sp 3 - 
hybridized bridging atom prevents the mixing of n- and cr- 
orbitals. This allows us to construct diabatic models using the 
ASD-CASSCF method. 

Finally the model Hamiltonians for hole and electron trans¬ 
fer processes of anti-[ 2.2](l,4)pentacenophane (PP) were cal¬ 
culated using ASD-RASSCF with N - 8. PP has re¬ 
cently been synthesized as a novel pentacene derivative.^The 
charge (especially hole) mobilities of large acenes are of in¬ 
terest for photovoltaic applications: 59 We used a molecular 
structure optimized by density functional theory (B3LYP/6- 
31G*)! 25 ®The S, T, and CT subspaces were included in the 
ASD expansion, in which RAS(7,4,7)[2,2] (7,4, and 7 orbitals 
were assigned to the RAS I, II, III spaces for each monomer, 
respectively; up to 2 holes and 2 particles were allowed in the 
RAS I and III spaces) was used for monomer states. We fixed 
the monomer Cl coefficients after five iterations and recom¬ 
puted the ASD-CI energies using the optimized orbitals (see 
Sec. IIA i. This procedure was iterated until consistency be¬ 
tween orbitals and monomer Cl coefficients was attained. The 
diabatic couplings obtained from ASD-RASSCF are 125.1 


and 11.6 meV for hole and electron transfer processes, respec¬ 
tively. 

IV. CONCLUSIONS 


We derived and implemented an orbital optimization 
method in the framework of ASD, in which a dimer wave 
function is expanded in products of monomer functions. The 
active-active rotations between ASD subspaces are included 
in the optimization. This not only removes the ambiguity in 
the choice of ASD subspace orbitals but also extends the ap¬ 
plicability of ASD to covalently linked dimers. The energies 
computed using both ASD-CASSCF and ASD-RASSCF con¬ 
verge rapidly with respect to the number of states in the ASD 
expansion. Since orbital optimization preserves the diabatic 
nature of the dimer basis functions, we were able to compute 
orbital-optimized model Hamiltonians for electron, hole, and 
triplet transfer processes. Work toward incorporating dynami¬ 
cal correlation on the basis of the ASD-CASSCF and ASD- 
RASSCF methods is in progress. The ASD-CASSCF and 
ASD-RASSCF methods also provide smooth potential energy 
surfaces and can be used to study reaction dynamics. All the 
programs developed in this work are available in the open- 
source bagel package.®! 


ACKNOWLEDGMENTS 


This work has been supported by Office of Basic Energy 
Sciences, U.S. Department of Energy (Grant No. DE-FG02- 
13ER16398). 




7 


1 H. van Amerongen, L. Valkunas, and R. van Grondelle, Photo¬ 
synthetic excitons (World Scientific, Singapore, 2000). 

2 G. R. Desiraju, Chem. Rev. 135, 9952 (2013). 

3 A. Nitzan, Chemical Dynamics in Condensed Phases (Oxford 
University Press, Oxford, 2006). 

4 S. Difley and T. Van Voorhis, J. Chem. Theory Comput. 7, 594 

( 2011 ). 

5 B. Kaduk, T. Kowalczyk, and T. Van Voorhis, Chem. Rev. 112, 
321 (2012). 

6 J. E. Subotnik, J. Vura-Weis, A. J. Sodt, and M. A. Ratner, J. 
Phys. Chem. A 114, 8665 (2010). 

7 E. C. Alguire, S. Fatehi, Y. Shao, and J. E. Subotnik, J. Phys. 
Chem. A118, 11891 (2014). 

8 B. R. Landry and J. E. Subotnik, J. Chem. Theory Comput. 10. 
4253 (2014). 

9 R. F. Fink, J. Pfister, A. Schneider, H. Zhao, and B. Engles, Chem. 
Phys. 343, 353 (2008). 

10 R. W. A. Havenith, H. D. de Gier, and R. Broer, Mol. Phys. 110. 
2445 (2012). 

11 W. Wu, P. Su, S. Shaik, and P. C. Hiberty, Chem. Rev. Ill, 7557 
( 2011 ). 

12 S. M. Parker, T. Seideman, M. A. Ratner, and T. Shiozaki, J. 
Chem. Phys. 139, 021108 (2013). 

13 S. M. Parker and T. Shiozaki, J. Chem. Theory Comput. 10, 3738 
(2014). 

14 S. M. Parker, T. Seideman, M. A. Ratner, and T. Shiozaki, J. Phys. 
Chem. C 118. 12700 (2014). 

15 S. M. Parker and T. Shiozaki, J. Chem. Phys. 141, 211102 (2014). 


E. Alguire, J. E. Subotnik, and N. Damrauer, J. Phys. Chem. A 
119, 299 (2015). 

G. K.-L. Chan and S. Sharma, Ann. Rev. Phys. Chem. 62, 465 

( 2011 ). 

G. Chaban, M. W. Schmidt, and M. S. Gordon, Theor. Chem. 
Acc. 97, 88 (1997). 

J. Olsen, B. O. Roos, P. Jprgensen, and H. J. A. Jensen, J. Chem. 
Phys. 89, 2185 (1988). 

T. H. Dunning, J. Chem. Phys. 90, 1007 (1989). 

See Supporting Information for the convergence behavior of 
ASD-CASCI and ASD-CASSCF with aug-cc-pVDZ basis set. 

G. L. Closs, P. Piotrowiak, J. M. Maclnnis, and G. R. Fleming, J. 
Am. Chem. Soc. 110, 2652 (1988). 

G. L. Closs, M. D. Johnson, J. R. Miller, and P. Piotrowiak, J. 
Am. Chem. Soc. Ill, 3751 (1989). 

R. Nogita, K. Matohara, M. Yamaji, T. Oda, Y. Sakamoto, T. Ku- 
magai, C. Lim, M. Yasutake, T. Shimo, C. W. Jefford, , and 
T. Shinmyozu, J. Am. Chem. Soc. 126, 13732 (2004). 

A. D. Becke, J. Chem. Phys. 98, 5648 (1993). 

P. C. Hariharan and J. A. Pople, Theor. Chim. Acta 28, 213 (1973). 

F. Weigend and R. Ahlrichs, Phys. Chem. Chem. Phys. 7, 3297 
(2005). 

R. Bula, M. Fingerle, A. Ruff, B. Speiser, C. Maichle-Mossmer, 
and H. F. Bettinger, Angew. Chem. Int. Ed. 52, 11647 (2013). 

J. E. Anthony, Angew. Chem. Int. Ed. 47, 452 (2008). 
bagel. Brilliantly Advanced General Electronic-structure Library. 
http://www.nubakery.org under the GNU General Public Li¬ 
cense. 


16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 



